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(57) Abstract: Improved systems and 
methods for coded aperture imaging of 
radiation-emitting sources. According to 
one aspect, the present invention is directed 
to reducing and/or eliminating the artifacts 
that are inherent in previous near field 
coded aperture imaging systems such that 
the improved sensibility and resolution of 
these systems can be practically utilized. A 
system and method of the present invention 
utilizes two projections of radiation from 
an object, the first through a first coded 
aperture mask pattern, and the second 
through a second coded aperture mask 
pattern, where a decoding array associated 
with the second mask pattern is the negative 
of a decoding array associated with the first 
mask pattern. Data from both projections is 
combined to produce a reconstructed object 
image that is substantially free of near-field 
artifacts. The present invention additionally 
relates to further improvements in the 
sensitivity and spatial resolution of coded 
aperture imaging applications. More 
particularly, the present invention relates to 
improvements in the design and fabrication 

of coded aperture masks for use in coded aperture imaging devices. Improvements can be made to the resolution, for instance, by 
selecting smaller pixel sizes for the opaque and transparent elements of the coded aperture mask. Also, it is possible to improve 
the signal-to-noise ratio for near-field applications by selecting an appropriate thickness for the mask. The signal-to-noise ratio 
may be further improved by selecting the appropriate mask pattern based on the particular characteristics of the near-field imaging 
problem at hand. 
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CODED APERTURE IMAGING 

BACKGROUND OF THE INVENTION 

Coded aperture imaging has been proposed in the past as a means for 
improving the spatial resolution, sensitivity, and signal-to-noise ratio (SNR) of 
5 images formed by x-ray or gamma ray radiation. For many imaging applications, 
coded aperture cameras have proven advantageous relative to other candidate 
systems, including the single pinhole camera and multihole collimator systems. In 
contrast to these other systems, for instance, the coded aperture camera is 
characterized by high sensitivity, while simultaneously achieving exceptional spatial 

1 0 resolution in the reconstructed image. 

In contrast to the single pinhole camera, coded aperture systems utilize 
multiple, specially-arranged pinholes to increase the overall photon transmission, 
and hence the sensitivity, of the imaging camera. In operation, radiation from the 
object to be imaged is projected through the coded aperture mask and onto a 

15 position-sensitive detector. The coded aperture mask contains a number of discrete, 
specially arranged elements that are either opaque or transparent to the incident 
photons. The raw signal from the detector does not reflect a directly recognizable 
image, but instead represents the signal from the object that has been modulated or 
encoded by the particular aperture pattern. This recorded signal can then be digitally 

20 or optically processed to extract a reconstructed image of the object. 

SUMMARY OF THE INVENTION 

Most commonly, coded aperture imaging techniques are employed in the 
field of astronomy. For instance, coded aperture techniques have been used for 
many years in the imaging of distant x- or gamma-ray astronomical sources. These 
25 applications are classic "far field" (i.e. small object angle) imaging applications, 
where, under the proper conditions, conventional deconvolution techniques maybe 
employed to produce ideal images of the source. 
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However, for near field applications, including nuclear medicine, molecular 
imaging, materials analysis, and contraband detection, among others, these same 
deconvolution methods do not produce ideal object images. Rather, in near field 
imaging applications, the reconstructed object images are corrupted by artifacts, 
5 including unsightly lines and/or bows through the object image. The presence of 
these artifacts render coded aperture systems far less attractive for near field imaging 
applications, particularly for applications such as nuclear medicine imaging, which 
rely heavily on visual inspection of the reconstructed object image. In the particular 
case of nuclear medicine imaging, for example, the presence of these artifacts may 

10 render it difficult or impossible to distinguish actual objects from artifacts in the 
reconstructed images, and may be a prohibitive factor weighing against the 
widespread adoption of coded aperture techniques. 

The present invention relates to improved systems and methods for near-field 
coded aperture imaging of radiation-emitting sources. According to one aspect, the 

1 5 present invention is directed to reducing and/or eliminating the artifacts that are 
inherent in previous near field coded aperture imaging systems such that the 
improved sensitivity and resolution of these systems can be practically utilized. 

In one aspect, the present invention is based upon the recognition that near 
field artifacts result primarily from the physical process of projecting radiation from 

20 the object through the coded aperture and onto the detector. It has heretofore been 
assumed that, for coded aperture arrays having perfect cross-correlation properties, 
the projected image obtained at the detector is the convolution of the object with the 
aperture mask. In fact, however, this projection contains a non-linear cos 3 (0) term 
affecting the image. For far field applications, including astronomical applications 

25 where the object-to-detector distance is many orders of magnitude greater than the 
detector size, this cos 3 (6) term can be approximated as equal to one, and the 
projected image on the detector is the straightforward convolution of the object and 
aperture mask. 

For the near field applications of the present invention, however, this cos 3 (0) 
30 term does not approximate as one, and therefore it gives rise to visible artifacts 
which corrupt the reconstructed object image. The present invention relates in one 
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aspect to systems and methods for reducing or removing these artifacts from the 
reconstructed image. More specifically, the invention comprises providing a 
radiation-emitting object within the field of view of a coded aperture camera; 
generating a signal corresponding to a first image of radiation from the object 

5 projected through a first coded aperture mask pattern; generating a signal 

corresponding to a second image of radiation from the object projected through a 
second coded aperture mask pattern, where the second pattern is the "negative" of 
the first; and then processing the combined data from these projections to obtain a 
reconstructed image of the object. 

10 In the context of the present invention, the second mask pattern is related to 

the first mask pattern because the second position is the one associated to the 
negative of the decoding pattern associated with the first physical mask pattern. For 
each coded aperture mask pattern, there is an associated decoding function, G, that is 
used to decode the recorded image and produce the source reconstruction. One may 

1 5 change the sign of G to produce -G, i.e. the negative of the decoding array. The 
second projection is then taken with the physical mask associated with -G, which 
according to this invention is referred to as the "negative" mask. As near-field 
artifact effects are dependent on G and -G, when the data sets from both projections 
are combined, the near-field artifacts cancel from the reconstructed image. 

20 Some mask families, including the URA and m-sequence families, have the 

particular property that the "negative" mask may be produced by simply 
interchanging the position of the opaque and transparent elements of the physical 
mask. This is because, for these families, the array, A, and decoding array, G, are 
identical. Thus, when signs change for the transparent and opaque elements of the 

25 decoding array, the elements simply interchange with one another. Also, because 
the array and decoding array are identical, the elements similarly change for the 
associated "negative" physical mask. The resultant mask, with the elements 
interchanged, is sometimes referred to as the "antimask" of the first mask. 

For some near-field applications, two physical masks maybe used to provide 

30 the two masks: the first being the mask pattern, and the second being the negative 
mask pattern. Depending on the particular mask family chosen, the second mask 
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could be the antimask of the first. Ia this example, the physical masks can be 
interchanged between projections, either by operator intervention, or by an 
automated process. 

In other cases, depending on the particular mask pattern used, a single 
5 physical mask may comprise both the mask and its negative mask. For some mask 
families, for instance, the negative of the original mask pattern (i.e. the mask 
associated with -G) may be obtained by simply rotating the physical mask by a given 
angle. For certain antisymmetric mask patterns, including a No-Two-Holes- 
Touching MURA pattern, for instance, a first mask pattern will produce its negative 

10 pattern when rotated by 90°. 

The use of a mask and its antimask has been reported previously in 
connection with far-field astronomical applications for the much different goal of 
reducing non-uniform background and compensating for systematic irregularities. 
For instance, in the context of a balloon or satellite-borne gamma-ray telescope, the 

15 mask-antimask technique may be useful to "average out' ' irregularities that may 
develop in the detector apparatus, particularly given the practical difficulties in 
accessing and actively maintaining the equipment. 

On the other hand, near field coded aperture imaging applications uniquely 
suffer from the problem of near field artifacts, a problem which may be minimized 

20 through the use of a mask and its negative, which may, in certain instances, include 
the use of the mask/antimask pair. Generally, it has been found that when the data 
sets of an image taken through a mask and its negative mask are combined, the 
object images reinforce, while the near-field artifacts cancel. Thus, the 
mask/negative mask technique may be easily and efficiently employed in near field 

25 problems to produce robust images that are substantially free from artifacts caused 
by near-field geometry. The present invention is particularly advantageous in that it 
is implemented primarily in hardware, it does not rely on heavy computation and it 
does not require any significant extension of the time needed to produce an image. 
In the context of the present invention, near-field imaging comprises coded 

30 aperture imaging applications wherein the projected image obtained at the detector is 
not the straightforward convolution of the object with the aperture mask, but is 
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additionally affected by a cos 3 (8) term that contributes to the reconstructed image in 
the form of visible artifacts. For many, but not necessarily all, near field 
applications, the object is less than about 10 meters from the detector. For other 
* applications, including nuclear medicine imaging, where the goal is often to collect 

5 the data as close to the source as possible, the object-to-detector distance can be less 
than about 1 meter. 

The present invention additionally relates to further improvements in the 
sensitivity and spatial resolution of coded aperture imaging applications. More 
particularly, the present invention relates to improvements in the design and 

10 fabrication of coded aperture masks for use in coded aperture imaging devices. 

Improvements can be made to the resolution, for instance, by selecting smaller pixel 
sizes for the opaque and transparent elements of the coded aperture mask. Also, it is 
possible to improve the signal-to-noise ratio for near-field applications by selecting 
an appropriate thickness for the mask. The signal-to-noise ratio may be further 

1 5 improved by selecting the appropriate mask pattern based on the particular 
characteristics of the near-field imaging problem at hand. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a schematic view of a coded aperture imaging system according to 
the present invention; 
20 Fig. 2A is an illustration of a coded aperture mask; 

Fig. 2B is an illustration of a coded aperture mask that is the negative of the 
mask of Fig. 2 A; 

Fig. 3 is a process diagram illustrating coded aperture imaging as a 
projection process and an object reconstruction process; 
25 Fig. 4A is an illustration of a 79 x 79 MURA coded aperture mask pattern; 

Fig. 4B is an illustration of a 77 x 77 Product Array coded aperture mask 
pattern; 

Fig. 4C is an illustration of a 62 x 62 No-Two-Holes Touching MURA coded 
aperture mask; 
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Fig. 4D is an illustration of a 63 x 65 m-sequence coded aperture mask 
pattern; 

Fig. 5 is a schematic view of the geometry of coded aperture imaging; 
Fig. 6 is an illustration of first order artifacts in a near-field coded aperture 
5 imaging application; 

Fig. 7 A is an image of a thyroid phantom projected through a mask pattern; 
Fig. 7B is an image of a thyroid phantom projected through the negative 
mask pattern of Fig. 7 A; 

Fig. 7C is an image of a thyroid phantom representing the sum of Figs. 7A 

10 andTB; 

Fig. 7D is an image of a thyroid phantom representing the subtraction of 
Figs. 7A and 7B; 

Fig. 8 is a schematic illustration of a point source being imaged with a 
mosaicked and a non-mosaicked coded aperture mask pattern; 
15 Fig. 9 is a plot of mask thickness versus a figure of merit for various mask- 

to-detector distances; 

Fig. 1 OA is an illustration of an asymmetric 11x11 No-Two-Holes- 
Touching MURA coded aperture mask pattern; 

Fig. 1 0B is an illustration of an asymmetric 23 x 23 MURA coded aperture 
20 mask pattern; and 

Fig. 1 1 is a perspective side view of two point sources in a three-dimensional 
space projecting coded aperture mask patterns of different sizes onto a detector. 

The foregoing and other objects, features and advantages of the invention 
will be apparent from the following more particular description of preferred 
25 embodiments of the invention, as illustrated in the accompanying drawings in which 
like reference characters refer to the same parts throughout the different views. The 
drawings are not necessarily to scale, emphasis instead being placed upon 
illustrating the principles of the invention. 
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DETAILED DESCRIPTION OF THE INVENTION 

In Fig. 1, a system 10 for imaging an object according to the present 
invention is illustrated, the system comprising a coded aperture 2, a detector 3, and a 
data processor 4. The coded aperture 2 and detector 3 define an imaging camera 5 
5 having a predetermined field of view (FoV). In operation, an object 1 or a portion of 
an object to be imaged is positioned within the FoV of the imaging camera 5, where 
the camera is at a selected distance from the object. Alternatively, the object can 
remain stationary, and the camera can be positioned such that the object or portion 
of the object of interest is within the FoV of the camera. 
10 The object emits radiation 6, such as, but not limited to, x-ray and/or y-ray 

radiation, that is modulated by the coded aperture 2 and impinges upon the detector 
3. As illustrated in Fig. 2A, the coded aperture generally comprises open transparent 
regions 1 1 and closed regions 12 that are substantially opaque to radiation that is 
emitted by the object. 

15 If the object is extended, it can generally be treated as comprising multiple 

point sources, each of which emits radiation. Each of these point-sources casts a 
particular shift of the coded aperture pattern on the detector 3. Thus, many different 
shifts, corresponding to the different point sources comprising the radiation emitting 
object, are superimposed on the detector. The detector 3 provides detection signals 

20 corresponding to the energy and pattern of the emitted radiation, and the processor 4 
can subsequently characterize the object 1 based upon the detection signals. 
According to one aspect, the processor 4 can characterize the object by 
reconstructing a visible image 7 of the object. In this embodiment, the imaging 
system additionally comprises a display 8 for illustrating the reconstructed object 

25 image. 

The detector 3 generally comprises a position sensitive detector for recording 
the transmitted emission signals. A single detector or a line detector can be used to 
record the spatial distribution of the transmitted emission signals by moving through 
the entire shadow-casting area within a plane. Preferably, however, the detector 
30 comprises a two-dimensional detector array, where the detection plane elements 

correspond to either a defined region of a continuous detector, or individual detector 
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units spanning the entire area in which the coded aperture casts a shadow. 
Conventional detector arrays, such as an Anger camera, including a Siemens E-Cam 
medical gamma camera, may be employed. 

One suitable detector arrangement, for instance, is to employ a scintillating 

5 material in conjunction with'photomultipliers. Useful scintillators include sodium 
iodide, cesium iodide, gadolinium oxyorthosilicate (GSO), bismuth germanate 
(BSO), and lutethium oxyorthosilicate (LSO). The scintillators can be constructed as 
discrete elements or as a single large crystalline sheet. An array of photomultipliers 
is typically used in conjunction with the scintillator material to accurately record the 

1 0 incident radiation, including its energy and position. Other detectors could include 
noble gas ionization chambers or solid state materials such as high purity 
germanium, cadmium zinc telluride (CZT), or cadmium telluride. 

The data processor 5, which can be a programmable computer, receives 
signals from the detector 3 corresponding to the detected radiation signals, and is 

15 programmed to characterize the object, including reconstructing an image of the 
object, based upon the transmitted radiation and the configuration of the coded 
aperture. 

The apertures can be in any shape, including polygons, circles, rings, or a 
combination of these shapes. Generally, the aperture shape is the same as the cross- 

20 sectional shape of the detector elements. 

hi certain embodiments, it may be advantageous to obtain multiple exposures 
of the object through one or more coded aperture masks. As described in further 
detail below, the problem of near-field artifacts may be greatly reduced where the 
coded aperture camera 5 obtains at least two images using different coded aperture 

25 arrays. More specifically, a first mask pattern can be used to obtain a first recorded 
image, and a second mask pattern that is the "negative," (and, in some cases, the 
"antimask") of the first can be used to obtain a second recorded image. The second 
mask pattern is called the negative of the first because the decoding array, G 2 , 
associated with the second mask pattern, is the negative of the decoding array 

30 associated with the first mask pattern, Q v In other words, -G 2 = G v The coded 
aperture can then take two pictures of the object, one through the mask and the other 
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through its negative mask, and the data processor can add the two data sets to 
eliminate near-field artifacts that would otherwise corrupt the reconstructed image. 

In certain embodiments, a single physical mask can serve as both the mask 
and its negative mask. For instance, depending on the axis of symmetry of the mask 
5 pattern of the coded aperture, a single mask may comprise a first mask pattern, and 
then the second, or Negative" of the first pattern when it is rotated about its center 
by a certain angle. Some mask patterns with an anti-symmetric axis of symmetry, 
for example, can be used as both the mask, and, when rotated by 90°, the negative 
mask. 

10 Two views of a coded aperture mask having an antisymmetric axis of 

symmetry are illustrated in Figs. 2A and 2B. The mask of Fig. 2A is a 62 x 62 No- 
Two-Holes-Touching MURA mask having a particular pattern of open elements 1 1 
and opaque elements 12. The mask is No-Two-Holes-Touching where no open 
element is directly adjacent to any other open element, and between each row and 

1 5 column of the mask array, there is a row or column of consisting entirely of opaque 
elements 13. In the "negative" mask, the open elements 1 1 switch places with the 
opaque elements 12, and the opaque strips 13 (which have values of zero in the 
associated decoding array) maintain their former location. Because of the 
antisymmetric axis of symmetry, this negative mask may be produced by simply 

20 rotating the mask of Fig. 2A by 90\ The result is the mask of Fig. 2B, which is the 
mask corresponding to the negative of decoding array of the mask of Fig. 2 A. 

The anti-symmetric mask is particularly efficient, as it does not require the 
fabrication of an entirely new mask in order to obtain the negative mask. Also, the 
use of an anti-symmetric mask simplifies the imaging process, as the only change 

25 needed to take the two images is to rotate the mask by 90°. This rotation process can 
be easily automated. 

However, there may be cases where two separate masks (i.e. mask and 
negative mask) are advantageous. For one, not every mask pattern has the desired 
axis of symmetry to enable the fabrication of an anti-symmetric mask. In cases 

30 where two separate masks are employed, the masks will need to be changed between 
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the two exposures. This can be done directly by the operator, or it could be done via 
an automated process. 

In Fig. 3, coded aperture imaging is illustrated generally as comprising a 
projection process and an object image reconstruction process. The projection 

5 comprises the physical steps of projecting the object through the aperture, to obtain a 
recorded image on the detector, as shown in step 14. The resulting recorded image, 
P, represents the correlation between the distribution of the activity on the object 
plane, O, and the aperture function at the mask plane, A, i.e. P = O x A. 

The reconstruction process generally entails decoding the recorded image 

10 computationally to produce an object reconstruction image, 6, as shown in step 15. 
The reconstructed image, 0, is produced by taking a periodic correlation between 
the recorded image, P, and the decoding array, G, (which consists of a data 
processing function derived essentially from the aperture mask pattern, A), such that 
6 = P ® G. The decoding array function is associated with the selected mask 

1 5 pattern, A, and can be programmed and stored in a digital computer. 

The combined process can thus be summarized by the following relationship: 

6 = P®G = (Ox A)®G = 0 *(A®G) = 0 *PSF (1) 

where * is the convolution function, and PSF represents the point spread 
function (i.e. the image that would be obtained if the object, O, was a point). It will 

20 be observed that if the PSF is a 8 function, then 6 = O, and the reconstruction is 
perfect. Thus, if one chooses A and G such that A® G = 8, then an ideal result is 
actually achievable. Several known families of aperture arrays have been found to 
exhibit this ideal behavior. A number of these families are discussed in R. Accorsi, 
Ph.D. Thesis in Nuclear Engineering, Massachusetts Institute of Technology, May 

25 2001, (and references therein), the entire contents of which are incorporated herein 
by reference in the entirety. 

Representative embodiments from a number of these families are illustrated 
in Figs. 4A-4D, including a 79 x 79 MURA mask (Fig. 4A), a 77 x 77 Product Array 
mask (Fig. 4B), a 62 x 62 No-Two-Holes-Touching MURA mask (Fig. 4C), and an 
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63 x 65 M Sequence mask (Fig. 4D). It should be noted for each family 
characterized by the ideal point spread function behavior described above, all 
patterns within that family also exhibit this ideal behavior, regardless of their 
physical dimensions. 

5 In the case of far-field imaging, such as the imaging of high energy x- and y- 

ray sources in astronomy, the mathematics of coded aperture imaging demonstrate 
that apertures with perfect correlation properties provide ideal point spread 
functions. However, applying the same techniques to near field problems, it has 
been found that the above deconvolution technique used for far field problems 

1 0 produce images that are affected by artifacts even where magnification is not pushed 
to achieve high resolution. 

In order to overcome the problem of artifacts in near field imaging, origins of 
these artifacts must be carefully examined. In particular, it has been found that the 
phenomenon of artifacts in near field images results primarily from the process of 

15 projecting the near-field object onto the detector. 

Fig. 5 is an illustration of the basic geometry of coded aperture imaging. 
Shown is an object plane 21, the aperture plane 22, and the detector plane 23. The 
number of counts recorded at a given detector position, ^ , in the detector plane is 

given by: 

20 P(r t ) oc JJ 0(f 0 )A^r 0 - f 0 )jcos 3 (8)rf 2 F 0 (2) 

with 6 = atgifc - f 0 1 / z) , O the distribution of the activity on the object 

plane, A the aperture function of the mask plane, and all other symbols explained in 
Fig. 5. If we define: 



l = --r 0 (3); 
a 
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0'(r) = 0^-^rj (4); and 



A\f) = A 



f <* 
-r (5) 



then Equation 2 can be rewritten as: 



P(f5)ocjj0'(|)^(fj-r 0 )cos 3 



V 



d 2 4 (6) 



/J 



5 O' and A' are, respectively, a rescaled and reflected form of the object, and a 

rescaled version of the mask pattern: the scaling coefficients are the magnification 
coefficients relative to the object and mask plane. Note that Equation 6 cannot be 
written as a correlation due to the near field term cos 3 (6). Therefore, the effect of 
near field is not a modulation of the object by a constant factor cos 2 (6), as in far field 

10 problems for off-axis sources. 

The reconstructed image is obtained from P via correlation with the decoding 
pattern, G\ associated to the mask A: 

6 = P®G' (7) 

1 5 where prime indicates that G was scaled according to A' . 

It will be apparent that if the object is far from the detector so that r i « z, 

we can make the far-field approximation: 



cos 3 (6)*l ( 8 > 
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10 



and write P as the convolution: 
P = O* * A' (9) 

which is equivalent to the correlation of the object (rescaled but not 
reflected) with the rescaled mask. Using the principles of coded aperture imaging 
described above, for ideal pairs (A,G), Equation 7 becomes: 

6 = (O' x A') ® G' = 9t[0' *(Q' ® A')] = &(P> *S) = 3t(0') (10) 

where 9t is the reflection operator. The reconstructed object is the object 
itself apart from a rescaling constant. However, if the far-field approximation, 
cos 3 (0) ^ 1, does not hold, then Equation 9 does not assume any simplified form and 
we do not reach this point. 

In near field imaging problems, such as nuclear medicine imaging, it is often 
important to collect data as close to the source as possible. This maximizes the 
sensitivity of the imaging system, and helps to minimize the dose levels of 
radioactive material given to the patient. Therefore, r i is often comparable to z and 

the far-field approximation does not hold. 

To understand the problem mathematically, the near-field term in Equation 6 
can be expanded in Taylor series to the second order. The result is: 




cos 3 



atg 



= cos 3 (^(%- 



z 



z 



J 
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1- 



3/z 2 



1+ 



_ flr la\-$ 2 5/2z 2 <z-^ 2 



b 2b 



1+ 



(11) 



This expansion is more accurate the larger the margin by which the condition 



J&_ = ^<1 (12) 



n 



10 



15 



is true. When high resolution is sought, magnification tends to be high, so 
a/b is generally low. Another way of looking at the same condition is to recognize 
that r. is a variable spanning the object, and, because of high magnification, if the 
object fits in the field of view, it is typically much smaller than the detector, and 
f ( > f 0 . Thus, the Taylor approximation is generally accurate. Indeed, it has been 

found that stopping at the second order is sufficient to explain the observed artifacts. 

It is possible to predict the occurrences of artifacts resulting from the 
contributions zero, first, and second order terms of the near-field equation, Equation 
11. More importantly, practical and effective systems and methods for reducing 
observable near-field artifacts in the reconstructed object image have been discovered 
in accordance with the present invention. 

Observing Equation 1 1, if the expansion is arrested at the zero order, one has: 



COS 3 ( at S( 



r. 



)) = C0S 3 (a#( m » (13) 
z 



which substituted in Eq. 6 gives: 
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P(^) = co f(atgfr)jjO\l)AXr t -&d 2 Z (14) 

The near-field effect is thus reduced to a prefactor modulating the image. As 
this term depends only on the detector coordinate, f t , it is relatively straightforward 

5 to correct for zero order artifacts exactly. After dividing the projection data, P, by 
this prefactor, the projection is again a convolution and the far field conditions of 
Equation 9 again apply. 

However, even correcting for the zero order effect, the overall problem of 
near field artifacts is not eliminated. Thus, it becomes necessary to correct for first 
1 0 and/or second order artifacts. 

The first order term is not negligible. From Eq. 1 1, its expression is: 



COS 



atg 



- a 7 



3 , >, L 31 z 

-cosVg(^)) 

z Ir. 



(i+^r) 



(15 



) 

1 5 Substitution into Equation 6 leads to, after zero order correction: 

z +r<l 4 

where ° is the scalar product. 

The integrand is the center of mass of the object "cut" by A' and the 
convolution gives the result as a function of the shift of A\ Assuming that A' covers 
20 uniformly the field of view and that the object is also reasonably uniform, the result is 
not a strong function of shift, and gives a constant contribution. The result can be 
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approximated by using a constant vector, O 1 , which is the center of mass of O' . 
Substituting in Eq. 9 thus gives: 




^G\r i + r})d 2 ri 



(17) 



For near field applications, the fraction in the integral can be substituted with 
5 r f I z 1 to further reduce to: 



6 oc o f o|J^G(^ +fi>/ 2 ^ (is) 



The second factor on the right-hand side of Equation 1 8 is the position of the 
center of mass of the open positions of G, the decoding pattern, as a function of the 
decoding position. This quantity depends upon the particular form of G. Fig. 6 

10 illustrates an example of these artifacts for a MURA mask pattern. Generally, as 
shown in Fig. 6, these artifacts will cause a vertical and/or horizontal line dividing 
bright and dark areas of the reconstructed image 31. As these artifacts depend on G 
only, they can be remedied through "pattern centering," that is, shifting the decoding 
pattern, G, (and thus the mask pattern, A), to move these lines to the borders of the 

15 image. For the case of the MURA shown above, this may be accomplished by 
shifting the decoding and mask patterns by half a period in both the x- and y- 
directions. The result is a corrected image 32, as shown in Fig. 6. Similar pattern 
shifting techniques may be used to mitigate first order artifacts in other mask array 
families, such as an m-sequence array or a product array. 

20 Note in Fig. 6 that the corrected image 32 still contains an artifact resulting 

from a linear modulation of the intensity making one comer of the image appear 
brighter than the others. However, this artifact can be eliminated completely with 
little effort. Returning to Equation 1 8, the first term on the right hand side of the 
Equation, (J , is the Center of Mass of the object. Thus, if object was centered on the 
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field of view, O' would be zero, and the term would cancel. Practically speaking, 

one can eliminate these artifacts by taking a first, raw image to estimate 0' and then 

making the necessary corrections before taking a second picture. 

We have observed that if the object is centered on the field of view, the first 
5 order artifacts disappear, but second order artifacts, normally hidden by the stronger 

higher-order artifacts, become visible. As illustrated in Figs. 7A and 7B, these second 

order artifacts for MURA patterns most commonly appear as semi-circular bows 41 

appearing on either side of the reconstructed object image. 

The starting point for analyzing the second order artifacts is substituting the 
10 second order terms of Equation 1 1 into Equation 6. For these artifacts, we must 

consider two terms. The first is: 

^)^-rA^JJ° ,( ^ (19) 

Here, the integrand is the (second) moment of inertia of O' (cut by A') with 
respect to an axis perpendicular to the object plane and passing through its center. If 
15 the open positions of A' are uniformly distributed (which is the case for the patterns 
shown in Figs. 4A-4D), the shift of A' does not modify the result. Upon decoding, 
the convolution of G with a constant gives a constant, and therefore this term can be 
neglected. 

The other second order term gives: 

20 p &) * 2 * , 2 2 p/J J o' (Ik en - ° £^ 2 5 (20) 

where r t is a unit vector. The integrand is the moment of inertia of O 5 (cut by 
A'), with respect to the axis on the plane of the object perpendicular to r ± . Therefore, 
it is a function with no dependence on (assuming, as before, that the open 



WO 02/056055 



PCT/US01/42420 



-18- 

positions of A' are uniformly distributed), and with an angular dependence of period 
w. The integrand can be rewritten as p/(/^), where p is the open fraction of A and 

/ is the previously discussed moment of O. The contribution to the decoded image 
is obtained by substituting in Equation 7: 

5 6 00 If : A^ frL)<7 G + ^ (21) 

This contribution is stronger the less isotropic the object, (i.e. if the object 
were a line source, for example). For an isotropic object, however, I is a constant and 
can come out of the integral. In that case, the artifact will again depend on G only. 
Equation 21 can thus be rewritten as: 

10 6 00 P J If 7T%|5^« + ^ (22) 

As a first approximation, one can ignore the denominator and conclude that 
second order artifacts have the shape of the moment of inertia of G (with respect to an 
axis perpendicular to its plane and passing through its center) as a function of shift. 
Not surprisingly, at second order the second moments of both O and G appear, and 
15 unfortunately, neither vanishes. However, for isotropic objects, the artifact depends 
only on a known quantity (G) and can thus be predicted independently of the object 

As previously discussed, these near-field artifacts can be easily and effectively 
reduced by recording two images of the object, one with a mask and one with its 
negative mask, and then adding the two recorded images. This technique relies on the 
20 fact that the near-field artifacts described above are dependent primarily on the array- 
decoding array pair, and not necessarily on the object to be imaged. 

More specifically, given a array-decoding array pair (A', G') which the ideal 
correlation properties discussed above, the pair (1-A', -G') likewise possesses ideal 
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correlation properties (except for a flat pedestal that does not distort the image and 
can be subtracted). The second mask is the negative of the first where the decoding 
array associated with the second mask changes signs relative to the decoding array 
associated with the first pattern. The sign changes cancel in the correlation to 
5 reconstruct the image, (see Eq. 10), so that the object images do not change sign. 
However, the first and second order artifacts predicted by Equations 18 and 22, 
depend only on G\ Therefore, when an image of the same object is taken with the 
negative mask, the artifacts change sign while the object image does not. Adding the 
two images reinforces the object image, and at the same time cancels the artifacts. 

10 This artifact reduction aspect is shown in greater detail in Figs. 7A-7D, which 

illustrates the experimental results of two exposures for a thyroid phantom. The mask 
exposure is shown in Fig. 7 A, and the negative mask exposure is shown in Fig. 7B. 
When the two images are added (Fig. 7C), the artifacts cancel out and the signal is 
reinforced. If the two images are subtracted (Fig. 7D), the opposite is true, and the 

1 5 artifacts reinforce while the obj ect cancels . 

This artifact reduction method has proven effective at reducing not only 
second order artifacts, but also the first order artifacts previously discussed. In future 
applications, such as applications requiring a wider field of view, higher order 
artifacts may additionally be encountered. It is believed, however, that the 

20 mask/negative mask technique described herein will be equally effective in reducing 
these higher order artifacts. 

The mathematics of near-field artifacts, and the artifact reduction techniques 
discussed above are also discussed in Roberto, Accorsi & Richard C. Lanza, <c Near- 
Field Artifact Reduction in Planar Coded Aperture Imaging," Applied Optics-IP, 

25 Volume 40, Issue 26, 4697-4705, September 2001, the entire contents of which are 
incorporated herein by reference. 

In addition to the near-field artifact reduction techniques described above, still 
further improvements to the coded aperture imaging systems/methods of the present 
invention may be realized by optimizing the parameters of the coded aperture mask. 

30 The field of view (FoV), resolution (A), and signal-to-noise ratio (SNR) are all 
dependent upon the physical dimensions and characteristics of the given mask. By 
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selecting the appropriate parameters for the mask, the performance of the coded 
aperture camera can be improved considerably. 

For instance, the field of view (FoV) may be increased by adopting a 
"mosaicked" mask pattenwe. replicating the mask pattern multiple times around a 

5 basic pattern, so that, for all points in the FoV, a shifted cycle of the basic pattern falls 
completely onto the detector. 

.The FoV is important in coded aperture imaging because the ideal point spread 
function behavior of an array-decoding array pair holds only if, for a given position in 
the field of view, the detector sees the shadow of a full mask pattern. Otherwise, the 

10 encoding is incomplete, and the periodic correlation with G does not produce a 8-like 
response. This is illustrated in Fig. 8, where a point source, P, is shown centered on 
the axis of detector 51, and projecting a shadow of the full mask pattern 52 onto 
detector 51. A second point source, P\ that is not centered on the axis of the detector 
51 is also projecting a shadow of the mask 52. However, this shadow does not fall 

1 5 completely on the detector. In this illustration, the field of view of the camera is 

limited to one point, P, the only point that casts a complete shadow onto the detector. 

In this example, the FoV could be improved by moving the mask closer to the 
detector until the shadows of all points of interest fall completely onto the detector. 
This solution is not ideal, however, as moving the mask closer to the detector will 

20 decrease the magnification, m, and consequently resolution will be lost. 

When a mosaicked mask pattern is utilized, however, the FoV may be 
increased without an attendant loss in resolution. Alternatively, for the same FoV, the 
resolution can be increased relative to the non-mosaicked mask. The only tradeoff is 
that a mosaicking the pattern requires a slightly more complicated mask. 

25 An example of a mosaicked mask pattern 53 is also illustrated in Fig. 8. 

Here, the basic pattern is held at the center of the mask and the copies are arranged 
around it. hi this embodiment, the left half is put to its right, and its right side to its 
left. This achieves a cyclic repetition as well as the more natural 2x2 lexicographic 
arrangement. Generally, though, the mosaicking technique cannot be extended to n x 

30 n mosaicking where n > 2. This would not increase the field of view over the 2x2 
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design, and would only corrupt the otherwise good field of view for the 2 x 2 portion 
of the mask. 

Another important parameter of the mask pattern is the mask pixel size, mp, 
which can be selected to boost the resolution of the coded aperture imaging system. 
5 For applications where the FoV and detector size are restricted or fixed, the 
magnification factor, m, is limited. Thus, one cannot rely indefinitely on the 
magnification to achieve the desired resolution. However, the resolution can still be 
improved by reducing the mask pixel size, mp, because the resolution may be defined 
as: 



Thus, reducing the mask pixel size will reduce (and, under this definition, 
improve) the resolution. 

There is, however, a limit to how small the mask pixel size.can be made. Two 
15 - factors which should be considered in attempting to optimize the resolution by 

reducing the mask pixel size include the detector pixel size and the mask thickness. 
Of these two, the most stringent restriction is set by the detector pixel size, dp, which 
is dependent not on the mask, but on the detector. 



20 (npm), the better the resolution, where: 
FoV = npmX (24) 

The limit on the npm is given by the product of npd, the number of detector 
pixels, and a, a sampling parameter; We have shown that the intensity of a point 
source is reconstructed independently of its position for any a ^2. Since: 

25 npd = a • npm (25) 



10 




(23) 



Consider that, for a given field of view, the higher the number of mask pixels 
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then: 

FoV = ^-Z (26) 
a 

which shows that for a given resolution, higher values of a lead to a smaller 
field of view, or, for a given field of view, a worse resolution. Consequently, a is 
5 preferably set to 2. 

The maximum value of npm gives the maximum field-of-view-to-resolution 
ratio. This value can be calculated from Equation 25 and from knowledge of the 
possible values of npm, which is restricted by the family chosen for the mask array. 
Once this value is known, the resolution (A) becomes fixed for a given field-of-view. 
10 By definition of m and a, it is also known that: 



mmp = a-dp (27) 

where dp is the detector pixel size (which is known because the detector is 
given). Since a is also known, this equation can be coupled with Equation 23 to 
obtain a system of two equations in two unknowns that can be solved for mp and m: 



15 



m 



a -dp -npm ^ ^ 

FoV + 
1 

mp = 



npm 



1 



FoV a -dp 



(28a-b) 



Thus, if one sets field of view = 9 cm, dp = 2.4 mm and npm « 62, selecting a 
mask with a pixel size of 4 mm will result in a resolution X=9 mm. Yet, by simply 
reducing the mask pixel size down to the smallest permissible value (which, 
according to the above equations results in a mask pixel size of 1.1 mm) the 
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resolution is improved from 9 mm down to X=1.45 mm, where the other parameters 
are held constant. 

Note that, in addition to decreasing the mask pixel size (mp), the resolution 
could also be improved by increasing the number of mask pixels (npm). For instance, 

5 if npm is 74, (which is a possible value for a self-supporting MURA mask), and all 
other factors are equal* the resolution in the above example can be improved to 1.22 
mm. However, a 74 x 74 self-supporting MURA cannot be made antisymmetric, thus 
requiring the fabrication of two masks in order to utilize the mask/negative mask 
artifact reduction process described above. As this example illustrates, one important 

10 limitation on the technique of increasing the number of mask pixels to improve 
resolution is the availability (or, rather, the lack thereof) of suitable mask patterns 
possessing optimal properties for near-field coded aperture imaging applications. 

Further improvements to the efficacy of the near field coded aperture imaging 
of the present invention, and specifically improvements to the signal-to-noise ratio 

15 (SNR) of the coded aperture imaging system, can be realized by selecting an 

appropriate mask thickness. Generally, when designing a coded aperture mask, the 
mask must be made sufficiently thick in order to properly encode the projected 
radiation. For instance, if the imaging system is measuring y radiation emitted from 
the object of interest, the closed regions of the mask must actually be opaque to most 

20 y rays. Undoubtedly, there will be some transmission through the mask. However, as 
the amount of this transmission increases, the SNR decreases, as evident from the 
following expression: 



SNR = 



^4B + 2(l + t 2 )I t 



(29) 



25 



with t = mask transmission; Sj = counts due to the reconstructed point; N = 
number of open mask pixels; B = counts due to uniform (unencoded) background; and 
I x = object intensity (i.e. the sum of all S^. 
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At the same time, any mask with finite thickness will invariably block at least 
some oblique rays that should pass through a mask opening. Further, as the mask 
thickness increases, less and less y rays will actually pass though, particularly in the 
case of near-field imaging where a substantial amount of the emitted radiation 
5 impinges on the mask at a non-zero angle. This effect is observable in the 

phenomenon known as 'Vignetting," where, as mask thickness increases, features of 
the projected image (starting with the small holes around the periphery of the image) 
begin to disappear. 

These factors (i.e. excessive transmission through opaque portions and 
1 0 "vignetting") should be considered in selecting a suitable material for fabricating the 
mask. Ideally, the mask material is one having minimum thickness for a given 
attenuation. In the particular case of y ray imaging, this is the material with the 
maximum product, \i x p, where \i is the attenuation coefficient at 140 keV and p is 
the material's density. Suitable practical materials for mask fabrication therefore 
15 include uranium (p,p = 48.97 cm 1 ), platinum (38.4 cm" 1 ), gold (35.9 cm" 1 ), tungsten 
(30.5 cm" 1 ), and lead (22.96 cm* 1 ). 

The mask thickness can be selected to optimize the tradeoff between 
transmission and "vignetting" described above. An appropriate mask thickness may 
be selected by first defining a Figure of Merit (FoM), which in this case may be the 
20 inverse of the standard deviation of the reconstruction of a point source (so that a high 
value is desirable) at the center of the field of view (so that pixel border effects are 
consistently handled), and, through simulation, graphing the FoM as a function of 
mask thickness. 

Another suitable Figure of Merit for the simulation is the signal-to-noise ratio 
25 (SNR) with which a point source at a given position within the field of view is 
reconstructed. This FoM can also be plotted as a function of mask thickness, as 
before. 

The object-to-detector distance, z, does have an effect on near-field artifacts, 
vignetting, and counting statistics, and should also be considered when selecting a 
30 suitable mask thickness for near-field coded aperture imaging applications. A graph 
of the FoM versus mask thickness for four values of the object-to-detector distance, z, 
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is illustrated in Fig. 9. Choosing a low value for z improves the counting statistics, 
but at the expense of increased near-field artifact and vignetting effects. However, by 
making z too large, there may be insufficient counting statistics for an accurate image 
reconstruction. As shown in Fig. 9, 30, 40, 50, and 60 cm values for z, all of which 

5 may be employed for near-field imaging applications, were used for the simulation. 

As illustrated in Fig. 9, for any given curve, the FoM initially increases with 
the mask thickness because the primary effect of increasing the thickness is to make 
the mask more opaque. However, as the thickness increases, vignetting effects begin 
to prevail, and, consequently, the FoM begins to decrease. The optimal mask 

10 thickness can be determined by finding the maximum value for a given object-to- 
detector distance, z. In this case, for z = 40 cm, for instance, the maximum FoM is at 
a mask thickness of approximately L5 mm. This value is not much different from the 
predicted maximum FoM values for the other object-to-detector distances, particularly 
the predicted FoM value when z = 30 cm. 

1 5 Tungsten is a particularly well-suited material for the coded aperture mask due 

to the attenuation characteristics noted above, as well as its reasonable cost, 
availability, and relative ease of fabrication. A tungsten mask with a thickness of 1 .5 
mm provides 99% attenuation (1% transmission) for y radiation at 140 keV. 

As previously discussed, there are a number of coded aperture mask families 

20 which exhibit the ideal imaging properties (i.e. A<8>G=6) for coded aperture imaging 
applications. Within these families, there are patterns having relatively little 
difference with respect to the achievable resolution and field of view. Also, the 
techniques for improved coded aperture imaging described above (i.e. mask/antimask, 
mosaicked pattern, reduction in pixel size, optimization of mask thickness, etc.), may 

25 be extended to any of these mask families. 

The choice of the mask pattern does impact one important characteristic of the 
imaging system: the signal-to-noise (SNR) ratio. One mask parameter in particular, 
the open fraction of a mask, p, (defined as the ratio of the open pixels to the total 
number of pixels) has previously been recognized as an important factor for SNR. In 

30 the case of a point source, for instance, the SNR is proportional to the square root of 



WO 02/056055 



PCT/US01/42420 



-26- 

the number of mask pixels, npm. Thus, it has generally been understood that the 
higher the open fraction of the mask, the better the SNR. 

It is also known, though, that the SNR advantage of a high open fraction 
becomes less pronounced as the object becomes less like a point source. In other 
5 words, as title object "concentration,"^ . , (defined as the ratio of the signal 

corresponding to the reconstructed point, S i5 to the total object intensity, I T ) decreases, 
the coded aperture overlaps data from the different radiation-emitting points. In the 
process of "undoing" this overlap between different sources, some SNR is lost. 

In addition to the object concentration, \|/ , , the number of "overlaps" further 

10 depends on p, the open mask fraction. The higher the open mask fraction, the more 
superimpositions are likely to occur. Consequently, there is a tradeoff in selecting a 
suitable open mask fraction: choosing a high value of p generally results in a high 
signal and thus improved SNR, but this also increases the superimposition of points 
on the detector, which undercuts or even eliminates this SNR advantage. 

1 5 There is thus a theoretical optimal open fraction for maximizing the SNR for a 

given coded aperture imaging application. However, as there are a finite number of 
available mask patterns to choose from, each associated with a certain value of p, it is 
unlikely that the theoretical optimal open fraction can be actually realized for a given 
coded aperture imaging application. Nevertheless, studying the effect of the open 

20 mask fraction on SNR can aid in the process of selecting a suitable pattern with 
sufficiently high SNR for a particular coded aperture imaging problem. One mask 
family, known as Uniformly Redundant Arrays (URAs), have previously been studied 
with respect to optimizing the open fraction to maximize SNR. Based upon these 
studies, it has been generally thought that relatively low values of p, less than 0.1, 

25 provide the best realistically achievable SNR. 

One aspect of the present invention is based upon the realization that the 
assumptions underlying the relationship between p and SNR for URA masks do no 
not have general application to all mask families. On the contrary, these assumptions 
regarding the relationship between p and SNR apply only to a particular subset of 

30 masks for which the mask pattern is equivalent to the decoding array, such as URAs. 
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For other masks, empirical observations reveal that the relationship between p and 
SNR is not the same as the case of URAs. For instance, out of the four candidate 
arrays illustrated in Figs. 4A-4D, the 77 x 77 product array (Fig. 4D) has the lowest p, 
approximately 0.09, (which is close to the predicted optimal open fraction for URA 
5 patterns), yet exhibits far inferior SNR than the other three arrays with larger open 
mask fractions: the 79 x 79 MURA (p = 0.5) (Fig. 4A), the 62 x 62 No-Two-Holes- 
Touching MURA (p = 0.125) (Fig. 4C), and the 63 x 65 m sequence (p - 0.5) (Fig. 
4D). 

The SNR can be generalized for all families of coded aperture masks, resulting 
10 in the expression: 



SNR iJ = 



p - oZ vm Z + <*+ #Z g I 



where N =. the number of open holes, t = the mask transmission, \\f = the 

concentration parameter, A and G the mask and decoding array, respectively, and § = 
background radiation. The present invention relates in one aspect to an improved 

1 5 method for coded aperture imaging wherein the mask pattern is selected in accordance 
with the above relationship to possess optimal or near-optimal SNR for the particular 
coded aperture imaging application. 

One family of mask patterns possessing advantageous SNR characteristics is 
the MURA family. A suitable MURA mask may be fabricated by securing a foil of a 

20 high-attenuation material, such as lead, to a supporting plate of a low-attenuation 
material, such as aluminum. The aluminum plate is intended to serve only the 
mechanical function of maintaining the integrity of the mask pattern, in particular, 
maintaining the positions of isolated opaque pixels of the mask. The mask itself can 
be manufactured using conventional techniques, such as milling. The MURA pattern 

25 can ideally be made antisymmetric, so that the single mask can serve as both the mask 
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and antimask for the purpose of near-field artifact reduction as described above. One 
example of an antisymmetric MURA is illustrated in Fig. 10B. 

Also, a self-supporting MURA mask may be employed. A self-supporting 
MURA is obtained from the standard MURA pattern by inserting opaque columns 

5 and rows between all the columns and rows of the standard MURA pattern. Thus, no 
opaque pixel is connected to any other opaque pixel by its corner only (which can be 
difficult to achieve from a fabrication standpoint), and all open pixels are separated 
from adjacent open pixels by at least one opaque pixel (i.e. No Two Holes Touching). 
The self-supporting MURA is characterized by a lower open mask fraction 

10 (approximately 0.125 as opposed to 0.5 for the standard MURA), and is further 
characterized by a slight reduction in SNR compared to the standard MURA. 
However, there are certain advantages to the self-supporting mask which may 
outweigh this slight reduction in SNR. For one, Tungsten can be employed as the 
mask material. As discussed above, Tungsten permits the fabrication of a mask 

1 5 characterized by high-attenuation at minimal thickness. Tungsten is an impractical 
material for a standard MURA, because it is very difficult to fabricate a Tungsten 
mask using conventional milling techniques. On the other hand, with a self- 
supporting mask, the support plate becomes unnecessary, and alternative fabrication 
techniques, such as photo-etching, can be employed. Tungsten is amenable to photo- 

20 etching, a known technique which utilizes acid to cut away the open holes of the 
mask. As with the conventional MURA, the self-supporting MURA can ideally be 
fabricated as an anti-symmetric mask, as illustrated in Fig. 10A. 

Other suitable materials for fabricating the mask include sintered Tungsten- 
based compositions (composed of greater than 90% Tungsten, for example). These 

25 materials have the additional advantage that they are more easily machinable than 
pure Tungsten. Suitable compositions for manufacturing coded aperture masks are 
commercially available under the trade names Densalloy™ and Heavimet™'. 

Once the mask pattern family is selected, the final parameter to consider is the 
size of the mask. This parameter is largely determined by factors previously 

3 0 discussed. For instance, assuming that the mask pattern is a self-supporting, No Two 
Holes Touching, mosaicked, pattern centered, antisymmetric MURA, the mask size is 
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determined by the detector parameters. From Equation 24, the number of mask pixels 
should be as high as possible to optimize resolution and field of view. However, the 
number of mask pixels is limited by the detector array characteristics, as well as the 
sampling parameter, a, which is set to two. 
5 If, for example, the detector array is 161 x 222 pixels, the maximum number 

of mask pixels is 161 x 16i; because all MURAs are square. However* since a = 2, 
this number is really only 80 x 80. Furthermore, since for self-supporting patterns 
every' other column and row is empty, the underlying MURA pattern can be, at most, 
40 x 40. Another characteristic of MURAs is that, for every pixel array n x n, n must 

10 be prime. Thus, the first available prime number that is less than 40 is 37, but a 37 x 
37 mask cannot be made antisymmetric. The next highest prime number n is 3 1, 
which will give an antisymmetric mask. Using an underlying 31x31 antisymmetric 
MURA pattern, and factoring in the opaque rows and columns for the self-supporting 
mask, the final pattern is therefore 62 x 62. 

1 5 It should be noted that a 37 x 37 MURA pattern may be employed (thus 

resulting in a 74 x 74 self-supporting MURA mask). This mask would exhibit 
improved resolution over the 62 x 62 mask due to the higher number of pixels. 
However, because this mask cannot be made antisymmetric, two masks (i.e. mask and 
antimask) would be needed in order to utilize the artifact reduction method of the 

20 present invention. 

The coded aperture systems and methods of the present invention may be 
advantageously used for high-resolution high-sensitivity imaging in nuclear medicine 
applications. Testing has been done with 1.15 mm (ID) capillary tubes injected with 
~ 1 50|iCi of 99m Tc, as well as 2d thyroid phantoms injected with ~200jiCi 99m Tc with 

25 the phantom at a distance of 40 cm from the detector. In the thyroid study, the two 
exposures (i.e. mask/negative mask) took approximately 8 minutes each. The results 
were in agreement with prediction as regards to near-field artifact reduction. The 
measured resolution was approximately 1.5 mm. 

In vivo studies have also been performed on mice injected with ~ lOOjiCi and 

30 -300jxCi of a bone agent, 99m Tc MDP, and similar tests have been done with a blood 
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agent, 99m Tc MPEG-PL. In these cases, high-resolution artifact-free images were 
obtained with improved SNR relative to the available technology, including 
collimators. In these studies, it has further been confirmed that the largest SNR 
advantages are realized when the object concentration, \j/ , is relatively large. The 

5 SNR can be increased, for instance, when a lead block is inserted to shield selected 
portions of the field of view, thus producing a high \j/ object of interest, and a much- 
improved reconstructed image. The present coded aperture imaging systems and 
methods maybe particularly beneficial for imaging small, "useftd" fields of view. 
Satisfactory results have also been achieved using higher energy (i.e> 140 
10 keV) isotopes. Although, in the case of a high-energy (511 keV) 18 F isotope, 

. achieving the same SNR with previously-described 62 x 62 MURA mask requires an 
exposure time approximately 50 times longer than the case of the lower-energy (140 

keV) 99m Tc isotope. 

New mask designs may be required to optimize performance for higher energy 

1 5 isotopes. For high-energy isotopes, a thicker mask with larger holes can be employed 
to achieve SNR performance equivalent to the lower-energy isotope examples using 
the previously-described mask, but at the cost of a worse resolution (e.g. 1.11 mm to 3 
mm). It is expected that the present invention can be extended for use with medium 
energy (e.g. m In, 67 Ga), as well as high energy (e.g. 131 I, 18 F, other PETs) isotopes. 

20 It is further expected that the coded aperture imaging according to the 

principles of the present invention will also prove advantageous for three-dimensional 
imaging applications, such as computer-aided tomography, or single photon emission 
computed tomography (SPECT). 

The coded aperture imaging system of the present invention can be used to 

25 obtain information about the depth of the imaged structures using a technique that is 
analogous to conventional x-ray laminography. According to this aspect, the coded 
aperture camera is able to image individual two-dimensional planes, or "slices", of the 
three-dimensional object. By consecutively imaging the individual slices making up 
the object, information regarding the object in all three dimensions may be obtained. 
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Classical x-ray laminography techniques are based on the relative motion of 
the x-ray source, the object to be imaged, and the detector. The x-ray source and the 
detector are moved synchronously (translationally or rotationally) in opposite 
directions. Because of this correlated motion, the location of the projected images of 
5 points within the object also moves. Only points from a particular slice of the object, 
the focal slice, will be projected always at the same location on the detector. Object 
structures above and below the focal slice will project at constantly changing points 
due to the relative motion of the source, object and detector. Consequently, the focal 
slice will be imaged sharply, while points outside of this slice will be superimposed as 
10 background intensity to the focal slice. Additional slices of the object can be imaged 
by displacing the object (either towards or away from the detector), and repeating the 
imaging process. 

In the context of the present invention, a method for laminograpy does not rely 
on any motion of the object, aperture mask or detector. Instead, the present invention 

15 utilizes the fact that point sources at different distances from the mask project mask 
patterns of different sizes onto the detector, where the size of the projected mask 
pattern is a unique function of the distance of the point from the detector. This is 
illustrated in Fig. 11, which shows two point sources in a three-dimensional space, 
each of which is a different distance from the mask and detector, projecting coded 

20 aperture mask patterns of different sizes onto the detector. After acquiring data from 
a three-dimensional object, one can scale the decoding pattern to match any size, and 
thus decode the data corresponding to any depth in the object. The decoding process 
will produce a sharp and artifact-free image of only the plane at the distance for which 
the decoding pattern is scaled. Point sources on other planes will be decoded 

25 incorrectly, and should blur into the uniform background. 

Moreover, as one can decode the same data after scaling the decoding pattern 
to another size, one can produce, from a single data set, images from different depths 
of the object. No object motion is required. 

Alternatively, the decoding pattern can be scaled to decode at a single depth 

30 only, and the object can be moved closer and farther from the detector to decode 
different focal slices of the three-dimensional object. 
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In addition to the nuclear medicine applications described above, the coded 
aperture imaging of the present invention may also be employed for the detection and 
imaging of radiation resulting from nuclear interrogation of a target object. For 
instance, coded aperture imaging maybe useful for the detection of contraband (e.g. 

5 explosives, drugs, and alcohol) concealed within cargo containers, suitcases, parcels, 
or other objects. Exemplary methods and apparatus for coded aperture imaging which 
employ nuclear interrogation of a target object are disclosed in United States Patent 
No. 6,205,195, "Coded Aperture Imaging," by Richard C. Lanza, the entire contents 
of which are incorporated herein by reference. 

10 In addition to nuclear medicine imaging and contraband detection, the 

principles of the present invention may prove useful for numerous additional coded 
aperture imaging applications, including materials analysis, scatter radiation imaging, 
and applications relating to detecting the movement or flow of radiation-emitting 
materials over time. 

1 5 While this invention has been particularly shown and described with 

references to preferred embodiments thereof, it will be understood by those skilled in 
the art that various changes in form and details may be made therein without departing 
from the scope of the invention encompassed by the appended claims. 
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CLAIMS 

What is claimed is: 

1 . A method of imaging radiation using coded apertures comprising: 

generating a first signal at a detector as a result of near-field radiation 
5 from an object being projected through a first coded aperture mask pattern; 

generating a second signal at the detector as a result of near-field 
radiation from the object being projected through a second coded aperture 
mask pattern, where the second mask pattern is associated with a decoding 
array that is the negative of a decoding array associated with a first pattern; 
10 and 

combining the first and second signals to obtain a reconstructed image 
of the object. 

2. The method of claim 1 wherein the object is located less than about 10 meters 
from the detector. 

15 3. The method of claim 1 wherein an image of the object projected through a 
coded aperture mask pattern is the convolution of the coded aperture mask 
pattern with the object that is modulated by a cos 3 (0) tenn, where cos 3 (0) is 
not approximately equal to one. 

4. The method of claim 1 wherein the emitted radiation is in the form of x-rays or 
20 gamma rays. 

5. The method of claim 1 wherein the emitted radiation is in the form of infrared 
radiation. 

6. The method of claim 1 further comprising interrogating the object with an 
energy source to produce the emitted radiation. 
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7. The method of claim 1 wherein the object is anatomical. 

8. The method of claim 7 further comprising evaluating the reconstructed image 
to determine the presence of tumors. 

9. The method of claim 7 wherein the obj ect contains a radioactive material. 

5 10. The method of claim 1 further comprising evaluating the reconstructed image 
to determine the existence of contraband articles. 

1 1 . The method of claim 1 further comprising: 

after radiation from the object is projected through the first coded 
aperture mask pattern, replacing a first coded aperture mask with a second 
10 coded aperture mask. 

12. The method of claim 1 wherein a single coded aperture mask comprises both 
the first mask pattern and the second mask pattern. 



1 3 . The method of claim 12 further comprising: 

after radiation from the object is projected through the first coded 
1 5 aperture mask pattern, rotating the coded aperture mask by a pre-determined 

angle to provide the second mask pattern. 

14. The method of claim 1 3 wherein the pre-determined angle is 90 degrees. 

15. The method of claim 1 wherein radiation from the object is projected through 
the first coded aperture mask pattern and radiation from the object is projected 

20 through the second coded aperture mask pattern for approximately equal time 

durations. 
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16. The method of claim 1 further comprising positioning the camera and/ or 
object such that the center of mass of the object is at about the center of the 
field of view of the camera. 

17. The method of claim 1 wherein the generated signals are processed using a 
decoding function associated with each mask pattern. 

1 8. The method of claim 1 7 wherein the decoding function is scaled to obtain a 
reconstructed image of a two-dimensional slice of a three dimensional object. 

19. The method of claim 1 wherein the generated signals from the first and second 
projections are added to cancel near-field artifacts that would otherwise appear 
in the reconstructed image. 

20. The method of claim 1 wherein processing data from the projections includes 
dividing the recorded data by a pre-factor to correct for zero-order artifacts. 

2 1 . The method of claim 1 further comprising: 

while obtaining an image projected through a mask pattern, changing 
the position of the object relative to the coded aperture mask and the detector 
to obtain an image of a cross-sectional slice of a three-dimensional object. 

22. The method of claim 1 wherein the mask patterns are mosaicked. 

23 . An apparatus for imaging an object emitting near-field radiation comprising: 

a detector array for detecting at least one of the energy or position of 
radiation emitted by the object and generating a signal in response thereto; 

at least one coded aperture mask, including a first mask pattern and a 
second mask pattern, disposed between the detector and the object, such that a 
first portion of the radiation is detected by the detector array after passage 
through the first mask pattern, and a second portion of the radiation is detected 



WO 02/056055 



PCT/US01/42420 



-36- 

by the detector array after passage through the second mask pattern, where the 
second mask pattern is associated with a decoding array that is the negative of 
a decoding array associated with a first pattern; and 

a data processor for processing signals, the signals generated in 
5 response to radiation passed through the first mask pattern and radiation 

passed through the second pattern, to produce a reconstructed image of the 
object that is substantially free of near-field artifacts. 

24. The apparatus of claim 23 wherein the detector array receives radiation from 
an object that is less than about 10 meters from the detector array. 

10 25. The apparatus of claim 23 further comprising display means for displaying a 
visible image of the object 

26. The apparatus of claim 23 wherein the emitted radiation is in the form of x- 
rays or gamma rays. 

27. The apparatus of claim 23 wherein the emitted radiation is in the form of 
15 infrared radiation. 

28. The apparatus of claim 23 wherein the at least one coded aperture mask 
comprises a first coded aperture mask having a first mask pattern and a second 
coded aperture mask having a second mask pattern. 

29. The apparatus of claim 28 additionally comprising means for interchanging the 
20 position of the first and second masks. 

30. The apparatus of claim 23 wherein a single coded aperture mask comprises the 
first mask pattern and the second mask pattern. 
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3 1 . The apparatus of claim 30 additionally comprising means for rotating a coded 
aperture mask having a first mask pattern by a pre-determined angle to provide 
the second mask pattern. 

32. The apparatus of claim 23 further comprising means for changing the position 
5 of the object relative to the coded aperture mask and the detector to obtain an 

image of a cross-sectional slice of a three-dimensional object. 

33. The apparatus of claim 23 wherein the mask patterns are mosaicked. 

34. A method for fabricating a coded aperture mask for use in a coded aperture 
imaging device, the coded aperture imaging device including an array of 

10 detector elements and having a selected field-of-view, the method comprising: 

selecting a pixel size for elements of a coded aperture mask, such that 
the size of the pixels approximately matches the smallest available pixel size 
to maintain the field-of-view of the coded aperture imaging device; 
selecting a coded aperture mask pattern; and 

15 on an aperture plate, defining a plurality of opaque and transparent 

elements, the opaque elements being substantially opaque to radiation emitted 
by an object to be imaged and the transparent elements being substantially 
transparent to radiation emitted by an object to be imaged, where the plurality 
of opaque and transparent elements is defined according to the selected coded 

20 aperture mask pattern, and each element has a size that is equal to the selected 

pixel size. 

35. The method of claim 34 wherein the elements are less than 2 mm wide. 

36. The method of claim 34 wherein the elements are approximately 1 .5 mm wide. 



37. 

25 



A method for fabricating a coded aperture mask for use in a coded aperture 
imaging device, where the device includes a detector for detecting radiation 
emitted from objects located at a pre-determined distance from the detector 
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and a processor for reconstructing an image of the object, the method 
comprising: 

selecting a coded aperture mask pattern; 

defining a figure of merit; 
5 through simulation, graphing die figure of merit as a function of mask 

thickness for the pre-determined object-to-detector distance; 

based upon the graphing, selecting a mask thickness that 
approximately matches the maximum value of the figure of merit for the pre- 
determined object-to-detector distance; and 

10 on an aperture plate having the selected mask thickness, defining a 

plurality of opaque and transparent elements, the opaque elements being 
substantially opaque to radiation emitted by an object to be imaged and the 
transparent elements being substantially transparent to radiation emitted by an 
object to be imaged, where the plurality of opaque and/transparent elements is 

1 5 defined according to the selected coded aperture mask pattern. 

38. The method of claim 37 wherein the mask thickness is approximately 1 .5 mm. 

39. The method of claim 37 wherein the figure of merit is the inverse of the 
standard deviation of the reconstructed image of a point source of radiation 
projected through a selected mask pattern, the point source located at the 

20 center of the field of view of the mask pattern. 

40. The method of claim 37 wherein the figure of merit is the signal-to-noise ratio 
of the reconstructed image of a point source of radiation projected through a 
selected mask pattern, the located at a given position within the field of view 
of the mask pattern. 

25 41 . A method for fabricating a coded aperture mask for use in a coded aperture 
imaging device, comprising: 

selecting a material for fabricating the coded aperture mask, the 
material having a thickness; 
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calculating the signal-to-noise ratio (SNR) for at least one coded 

aperture mask pattern using the equation: 

(l-t)N\\f if 

SNRjj = — Yjz 



where t is an estimate of an amount of radiation from an object to be imaged 
5 that is transmitted through opaque elements of the mask, N is a number of 

open holes in the mask, \j/ is an estimate of the concentration of the object to 
be imaged, A is a mask array function, G is a decoding array function 
associated with the mask array function, and £ is an estimate of background 
radiation; 

1 0 based upon said calculation, selecting a coded aperture mask pattern; 

and 

on the selected material, defining a plurality of opaque and transparent 
elements, the opaque elements being substantially opaque to radiation emitted 
by an object to be imaged and the transparent elements being substantially 
15 transparent to radiation emitted by an object to be imaged, where the plurality 

of opaque and transparent elements is defined according to the selected coded 
aperture mask pattern. 
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